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Preface 


Stephan  Prokofievitch  Timoshenko  was  bom  on  December  22,  1878  in  the 
Ukraine,  After  making  substantial  contributions  to  applied  mechanics  in 
Russia,  he  left  Russia  in  1920  for  Serbia,  and  in  1922  emigrated  to  this 
country.  After  five  years  in  industry,  (Vibration  Specialty  Company  and 
Westinghouse) , he  joined  the  faculty  of  the  University  of  Michigan  as  a 
professor  of  graduate  mechanics  and  later  (1936)  accepted  an  appointment 
at  Stanford  University. 

The  subject  of  this  study.  Traveling  Loads  on  the  Timoshenko  Beam, 
seems  particularly  appropriate  in  this,  the  centenary  of  Timosheko's 
birth.  In  addition  to  developing  the  improved  beam  theory,  (Philosophical 
Magazine,  1921-1922)  which  now  bears  his  name,  Timoshenko's  childhood 
interest  in  railway  engineering  was  fulfilled  through  many  works  on  the 
response  of  rails  and  bridges  to  moving  and  pulsating  loads.  Equation  110 
of  this  paper,  the  response  of  a Bernoulli-Euler  beam  to  a moving  force 
was  given  by  Timoshenko  in  the  Bulletin  of  the  Polytechnical  Institute  cf 
Kiev  in  1908  and  is  included  in  his  Vibration  Problems  in  Engineering. 

In  view  of  Professor  Timoshenko's  interest  in  this  subject  it  is 
somewhat  surprising  that  the  comparisons  made  here  were  not  given  long 
ago  by  Timoshenko  himself. 
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Abstract 


A transverse  force  travelling  along  an  Infinite  string  or  a beam  at 
critical  values  of  constant  velocity  generates  unbotuided  amplitudes.  In 
the  absence  of  dissipation.  This  resonance  Is  analagous  to  the  unbounded 
amplitudes  generated  by  a stationary  force  oscillating  at  one  of  the 
natural  frequencies.  The  response  of  a finite  elonentary  beam  to  a moving 
force  of  constant  amplitude  can  be  determined  In  terms  of  the  eigen- 
functions of  the  beam.  Modification  of  elementary  beam  theory  to  take 
into  account  the  effects  of  rotatory  Inertia  and  shear  leads  to  the 
Timoshenko  beam  theory,  from  which  a new  set  of  eigenvalues  and  eigen- 
functions can  be  determined.  These  eigenfunctions  can  be  shown  to  have 
an  orthogonality  relationship  which,  although  unusual,  permits  the  solu- 
tion of  Initial  value  and  non-homogeneous  problems.  The  procedure  for 
solving  such  problems  Is  given,  and  applied  to  the  problem  of  a travelling 
load  on  a finite  Timoshenko  beam  with  arbitrary  end  conditions.  Results 
are  obtained  for  the  case  of  pinned  ends,  and  compared  with  those  from 
elementary  theory.  Results  of  particular  significance  are  that  the  distri- 
bution of  critical  speeds  Is  altered  significantly  through  Inclusion  of 
rotatory  and  shear  effects,  and  that  a shear  wave,  not  present  In  the 
results  from  elesientary  theory,  can  be  Identified  and  is  shown  to  play  a 
major  role  In  determining  the  response. 
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I.  Introduction 
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The  Timoshenko  beam  has  been  considered  in  many  studies,  and  has 

found  numerous  applications  since  it  was  first  presented^.  Coiq>arisons 

have  been  made  between  the  predictions  of  this  more  exact  theory  and  the 

predictions  of  Bernoulli-Euler  theory  for  the  natural  frequencies  of  finite 

beams,  and  for  the  transient  response  of  beams  of  infinite  or  semi-infinite 
2 3 

length  ’ • While  a number  of  techniques  have  been  applied  to  the  solution 
of  problems  involving  the  Timoshenko  beam,  eigenfunction  expansions  do  not 
appear  to  have  received  much  attention  in  the  solution  of  initial  value 
problems,  and  in  the  determination  of  the  forced  response. 

It  is  evident  that  the  eigenfunctions  of  the  Timoshenko  beam  are 
orthogonal  in  the  usual  manner  in  the  case  of  simply  supported  ends.  Al- 
though a non-classic  orthogonality  for  other  boundary  conditions  has  been 
demonstrated^’^,  and  it  has  been  shown  that  arbitrary  initial  conditions 

5 

can  be  s.*tisfied  by  a series  of  the  eigenfunctions  , these  results  do  not 
appear  in  the  standard  texts  on  structural  vibrations.  It  also  does  not 
appear  to  be  well  known  that  the  response  to  a distributed  time  dependent 
force  can  be  written  in  terns  of  the  eigenfunctions.  Thus  it  appears  appro- 
priate to  set  down  a formulation  of  the  orthogonality  condition  idiich  leads 
directly  to  a procedure  whereby  initial-value  and  prescribed  force  prbblena 
for  the  Timoshenko  beam  may  be  treated.  This  procedure  will  then  be 
applied  to  two  example  problems,  the  transient  response  to  a uniform  im- 
pulse, and  to  the  subject  problem. 
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II.  Anal;>  is  of  Timoshenko  Beam 
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The  development  of  the  equations  of  motion  will  first  be  outlined, 
primarily  as  a means  of  defining  the  terms  to  be  employed.  The  deflections 
of  the  beam  are  taken  to  be 

u^(x,z,t)  - - z 4,  (x,t)  (1) 

u^(x,z,t)  - a)(x,t)  (2) 

plus  whatever  small  displacements  arise  through  Poisson  contractions.  Such 
displacements  will  be  presumed  to  give  rise  to  no  strain  energy,  nor  any 
kinetic  energy.  The  beam  geometry,  and  the  sign  conventions  for  displace- 
ments and  loads  are  depicted  In  Figure  1.  It  Is  to  be  noted  that  a 
distributed  moment,  not  present  In  elementary  theory  has  been  provided 
for,  although  It  Is  not  likely  to  occur  In  practice. 

From  the  figure,  and  Equations  1 and  2, 

i e^(x,z,t)  - - (3) 

i (m  (4) 

The  conmia  here  denotes  partial  differentiation  with  respect  to  the  coordi- 
nate. 


Applying  the  stress  strain  lav  leads  to  a uniform  shear  stress  which, 
when  Integrated  over  the  section^ leads  to  a pseudo  shear  force 


V 


rr 


J ] w(«J^3j-<(')dA 
A 


pA(w  -Ip) 


(5) 


with  u being  the  shear  modulus.  Since  the  actual  shear  stress  Is  not  uni- 
form, It  Is  customary  to  Introduce  a correction  factor  k,  leading  to  an 
expression  for  the  actual  shear  force, 

V(x,t)  ■ picA(u)  -♦)  (6) 
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Various  eatloates  of  the  shear  correction  factor  have  been  given  . 


Taking  the  energy  due  to  the  shear  to  be  the  product  of  the  average  shear 

stress,  V/A,  and  the  uniform  shear  strain,  the  total  strain  energy  in 
the  beam  is  found  to  be 

1 O t?  9 

(7) 


u - j^l  + |(z*^^)^}dAdx 


o A 

The  length  of  the  beam  Is  taken  as  L,  and  the  Young's  modulus  to  be  E. 

The  kinetic  energy  of  the  beam  Is 

L 


K ■ I I ^7  7 (**)^}dAdx 


(8) 


o A 

and  the  potential  due  to  distributed  forces,  moments  and  end  loads  Is 

L 

A - - I {qu  + mt|>}dAdx  - |m<|»  + Vul^  (9) 

o 

where  q amd  m are  the  distributed  force  and  moment,  per  unit  length, 
respectively.  M and  V are  resultant  moments  and  shear  forces,  applied  at  the 
beam  ends.  Performing  the  area  Integrals  and  substituting  Into  Hamilton's 
principle  yields 

.2 


.^2  r.L 


. f ,UkA  . . El  ^ 2 pA  .2  pi  .2 

' J J 'T-  2“ 


n r u 

L-  « qu  _ iBi|i}dx  - IM)()  + Vw]^  dt  = 0 

(10) 

where  1 ■ 

t 2 

j z dA.  Taking  the  variation.  Integrating  by  parts,  and 

collect- 

Ing  terms 

leaves  two  Euler  equations  and  two  boundary  conditions. 

" ®^’l',xx  ■ u<A(m^jj-1i)  + pi*  - m - 0 

(11) 

- yicA(«i)  ^-«i)  + pAw  - q - 0 

**  *x 

(12) 

!rutcA<(i»  - 0 

•*  o 

(13) 

|(EI*  - 0 

pA  O 

(14) 

Modulus,  area^and  moment  of  Inertia  have  been  taken  as  uniform,  although 
the  extension  to  nonuniform  sections  Is  not  difficult.  Thus,  satisfying 


i 


a prescribed  shear  force  or  displacenent  and  moment  or  rotation  at  each 
end,  together  with  satisfaction  of  the  differential  equations  on 
0 < X < L leads  to  an  extremum  of  the  time  integral  of  the  Lagrangian 
function. 


Ill,  Solution 


It  Is  convenient  to  put  Equations  11-12  In  a dlaenslonless  form. 
To  this  eno,  we  Introduce  a dimensionless  length, 


n ■ x/r 

where  r • I/A 

(15) 

a dlaenslonless  time 

T ■ tc/r 

where  c^  ■ yic/p 

(16) 

and  dlaenslonless  displacements,  moment  per  unit  length,  and  force  per 
unit  length  through  the  definitions 


♦ ■ w/r 
M • m/picA 
Q - qr/jjicA 

The  elastic  properties  can  then  be  described  by  a single  quantity 

^ KU  K 


(17) 

(18) 
(19) 


(20) 


The  result  of  these  substitutions  Into  Equations  11  and  12  are  dlaenslon- 
less forms  of  the  governing  differential  equations. 


7*  TT  ® 

* ni)  fX)  • TT 


-♦  +Q-0 

*n»1  tn  »TT 

The  homogeneous  boundary  conditions  become: 

(4  - - <*)  or  ^ must  vanish,  and 


tT) 


. or  p must  vanish 
#n 


(21) 

(22) 

(23) 

(24) 


at  n * 0 and  n = L/r  i.It  Is  convenient  to  rewrite  Equations  21  and  22 
as  a single  matrix  equation, 

(25) 


AU  +BU  +CU-DU  +P-0 
inn  ,n  itt 


where 


U 


(26) 
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F ■ { } and 

'q' 


y(5]  F® 

A=  .B=  .C=  .D=  (28) 

0 1 • 0*  00*  01  v^'^/ 

We  first  address  the  homogeneous  problem,  F = 0.  Let 


'f(n)e(T) 
♦ (n)0(T) 


u(n)e(T) 


Here  U Is  a vector  of  the  two  Independent  displacements,  but  U^U  is 
scalar  and  ve  may  obtain  a separation  of  variables  by  division. 


= - X 


The  differential  equation  to  be  satisfied  by  the  eigenfunctions  is 

AU  + BU  + CU  + ADU  - 0 (31) 

*nn  tb 


All  homogeneous  boundary  conditions  satisfying  equation  23  and  24  may  be 
put  into  the  form 

V,n  + 82^^  * ° 

and  must  be  satisfied  at  0 and  t.  The  necessary  matrices  B^  are: 

For  a clamped  end:  Fo  Ol  Fl  Ol 

•1  ■ [o  »J  • ^ L-  J 

For  a simply  supported  end: 

For  a free  end: 

"l  ol  To  0 

Bn  - , B - (35) 

1_0  ij  ^ [_-l  0_ 

and  for  a guided  end: 


( 


Satisfaction  of  the  matrix  differential  equation  (31)  with  a boundary 
condition  for  each  end  selected  from  33-36  permits  the  eigenfunctions  to 
be  determined. 

A.  Wave  Propagation 

The  propagation  o£  a sinusoidal  wave  train  in  a beam  of  infinite 
length  has  been  the  subject  of  a number  of  studies  . The  motion  is  then 
assumed  to  be 

ti  ■ {^>  exp  j (kx  - ot)  (37) 

B 

Ih  terms  of  the  dimensionless  variables  used  here,  we  have 


“ t^^expj  (krn  - — t) 
♦ B c 


(38) 


Satisfaction  of  equations  21  and  22,  with  M and  Q set  to  zero  requires 


that 


2 2 

yA  (-k^r^)  + Bjkr  - A + A -2-5-  - 0 


(39) 


B(-k^r^)  - Ajkr  + BH^r^/c^  - 0 (40) 

From  these,  the  dispersion  relationship  or  frequency  spectrum  is  obtained; 

0 (41) 


- kV)  - kV 

c* 


and,  if  desired,  the  ratio  A/B. 

In  terms  of  a dimensionless  frequency,  > ar/c^and  a dimensionless 
wavenumber,  K - rk,  this  becomes 

O A 

(42) 

For  any  K (wavenumber,  or  reciprocal  wavelength),  two  values  of  frequency 
are  found.  This  dispersion  relationship  is  shown  (for  y « 4)  in  Figure  2 
as  the  solid  lines.  For  the  lower  branch,  and  small  values  of  K,  it 
can  be  shown  that 


- (];  + (1  + Y)*i^Jl  + yK^  - 0 


X = yK 
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DIMENSIONLESS  FREQUENCY 


which,  in  terms  of  there  same  parameters,  is  the  frequency  spectrum  for 
the  Bemoulli-Euler  beam.  As  can  be  seen  in  Figure  2,  the  deviation 
between  the  spectrum  for  the  Timoshenko  beam  and  the  Bemoulli-Euler  beam  begins 
to  become  significant  at  K ^ 0.2,  In  the  Infinite  beam,  no  boundary 
conditions  need  be  satisfied,  and  motion  Is  possible  at  any  point  on  the 
frequency  spectrum.  We  will  see.  In  later  sections,  that  the  Introduction 
of  boundaries,  as  In  a finite  beam,  restricts  the  possible  isotlons  to 
discrete  points  on  this  same  frequency  spectrum.  i 

B,.  Orthogonality  of  the  Eigenfunctions. 

Let  us  denote  by  a set  of  functions,  each  of  which  satisfies  one 

of  the  allowed  boundary  conditions  at  n^O,  and  one  (not  necessarily  the 

same  one)  at  For  each  of  these  functions,  there  Is  an  eigenvalue 

X.,  which  we  assume  to  be  distinct, 
n 

Thus,  for  any  n,  ve  have 

+ BU“^  + C0“  + X DU“  • 0 (44) 

For  some  other  eigenvalue,  X^^  ^ X^,  we  have 

AU®  +BU“  +ai"+XDU“-0  (45) 

,nri  ,n  m 

Multiplying  the  first  by  and  the  second  by  where  T denotes  the 
transpose  of  the  2-vectors,  integrating  over  the  length,  and  subtracting, 
we  have 

1 1 
I - u“’^AU^„„}dn  + 1 {0“'^BU%  - u“’fB0^„}d!i 

+ I {0"’'‘co“  - u“’^C0"}dn+  I - X„u“’‘DU'“}dn  - 0 (46) 

o o 

where  1 * L/r,  Integrating  the  first  by  parts  leaves  only  boundary  terms, 
since  A Is  symsietrlc.  Integrating  only  the  second  t«m  of  the  second 


10 


Integral  by  pates  leaves  only  boundary  terms,  since  B la  skew  syoaetrlc 


Since  C and  D are  symmetric,  the  reductions  are  evident.  We  find: 

I 

lU®'^AU“  - U“'‘^An*  - - (X  -X  ) [u®’'DU“d  (47) 

,n  ,n  'o  m n j h 

o 

Evaluation  of  the  boundary  term  is  done  more  easily  (and  more  convincingly) 
by  returning  (through  Equations  28  and  29)  to  the  components  of  the  matrix 
U.  We  find  the  boundary  terms  to  be 


T 0 
0 1 


<» 

J .nj 


Y 0 
0 1 


f ,m 
'I'  _ 


1 


0 

-1 


.n  ,n  ,ri  ,n  'o 


(48) 


But,  each  product  vanishes  if  any  of  the  four  boundary  conditions  given 
by  Equations  13  and  14,  or  by  Equations  33-36,  are  satisfied  at  each  end. 
Thus,  for  X^  ^ X^,  we  have  the  desired  orthogonality  condition  in  terms 

of  the  vector  U t 

[ u“^DU"Jn  - 0 , (49) 


or,  in  terms  of  the  components. 


J + ♦“♦“}dn  ■ 0 , when  n m. 

In  terms  of  the  physical  rotation  and  displacement,  the  orthogonality 


(50) 


condition  may  be  rewritten 


{T“IT“  + W*AW“}dx  - 0 


vher« 


n n 

W « r f 


(51) 

(52) 


4.5 


An  equivalent  version  of  this  result  has  been  given  ' , but  does  not  seem 
to  be  widely  known. 


Although  only  honogeneoua  boundary  conditlona  have  been  dlacuaaad. 

It  la  eaally  ahown  that  orthogonality  exlata  (In  the  aaae  aenae)  If  lineal 
or  toralonal  springe  are  applied  at  the  enda*  but  does  not  1£  trsnslatory 
or  rotatory  masses  are  attached  to  the  end(s). 

Thus,  orthogonality  between  the  eigenfunctions  of  the  Timoshenko  beam 
exists  as  it  does  in  the  Bernoulli-Euler  beam,  if  the  eigenfunctions  are 
viewed  as  a 2-vectotsof  dlaplacement  and  rotation,  rather  than  as  dis- 
placements alone. 

C.  Initial  Value  Problems. 

We  now  asstime  that  a Timoshenko  bean,  having  appropriate  boundary 
conditions  at  each  end.  Is  released  at  time  zero  Into  free  vibration, 
with  specified  Initial  conditions 

♦(x,o)  - f(x) 

(i)(x,o)  - g(x) 

(53) 

i|>(x,o)  - h(x) 
u(x,o)  - k(x) 

We  look  for  a separable  solution.  In  the  form  of  Equation  29,  employing 
eigenfunctions,  given  by  44,  for  each  eigenvalue,  1^.  Solutions  to  the 
differential  equation  and  boundary  conditions  exist  only  at  discrete  val- 
ues of  the  wavenund>er,  R of  the  frequency  spectrum  for  the  Timoshenko  Beam, 
as  noted  In  Figure  2.  For  each  wavenumber,  however,  two  natural  frequencies 

are  fotind,  and  may  be  denoted  /X  , and  /X  ,,  respectively.  This  is 

ii  g X n g b 

demonstrated  by  the  dashed  lines  of  Figure  2.  Thus,  for  each  wavenumber 
k , m 3 1,2  •,  we  will  find  two  eigenvalues  X,  idiich  we  may 


number 


n 2m-l 


^ 4 


X = X„  = X - 
n 2m  m,2 

This  sequence,  while  convenient,  does  not  order  the  eigenvalues  in  ascending 
magnitude,  as  is  evident  from  Figure  2.  While  two  eigenvalues  arise 
from  the  same  wavenumber,  they  do  not  generate  the  same  eigen- 
function, for  the  ratio  ^ Consequently, 

as  will  be  more  evident  from  examples  to  be  given. 

We  let  n 

“ H'Cn)  \ 

tt  - ^ ^ 

n*l  •(!!) 


{ } satisfies  Equation  44.  The  satisfaction  of  initial  con- 


where  U 

dltlons  (Equation  53)  requires 


i|)(n,o)  “ f(n)  ■ I 
n“l 

09 

u)(ii,o)  - g(n)  - r^  a *”(n) 
n-1 

09 

4i(n,o)  - hCn)  ■ I v^Cn) 

n*l 

00 

u)(n,o)  ■ k(n)  ■ rl  — b •°(n) 
n-1  “ 


(56) 

(57) 

(58) 

(59) 


Multiplying  Equations  56  and  58  by  t".  Equations  57  and  59  by  $ , adding 
56  to  57  and  58  to  59,  integrating  over  the  beam  length,  and 


applying  the  orthogonality  condition.  Equation  50,  we  find 

J {f(n)’f"(n)  + ■“  g(n)9*(n)}<*n 

^ _ Jo g 


b. 


i [{(»“)'+  (•“)>dn 

I {rh^)'r®(n)  + k(n)*"<n)}dn 


® * /—  Z~2  ir~5 

cA^  {(T®)  + (♦“)  }dn 


(60) 


(61) 
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Thus,  the  necessary  coefficients  for  the  expansion  given  by  Equation  55 
for  the  initial  value  problcD  specified  by  Equations  53  have  been  deter- 
mined. An  equivalent  result  has  been  previously  given^. 

0.  Forced  Vibrations. 

We  now  return  to  the  general  problem,  Equation  25,  of  a beam  subjected 
to  an  arbitrary,  time  dependent,  load  and  moment  distribution.  Esch  end 
of  the  beam  is  assumed  to  satisfy  an  appropriate  homogeneous  boundary 
condition.  Only  a particular  solution  is  required,  as  the  method  of  the 
previous  section  may  be  used  to  furnish  the  homogeneous  solution. 

We  first  obtain  the  eigenfunctions,  and  expand  the  load  distri- 
bution and  the  dynamic  response  in  these  functions,  as 


and 


y - I D F„(T)0"(n) 
n-1 

m 

U - I G_(T)u®(n) 
n*l 


(62) 


(63) 


Substituting  into  Equation  25,  multiplying  by  and  integrating  over  the 

length,  we  find 
I 

\ tZ  + I tj“^BG^(T)U^^  + I o"’^CG^(T)u“}dn 

o j ^ “ (64) 

- I I 0"’^DG^'(T)u“dn  + jl  D*^^n(T)  0 u”dn  - 0 

o ° o'* 

The  first  three  terras  may  be  evaluated  by  recognizing  that  they  contain 

the  first  three  terns  of  Equation  44.  We  find 
t 

[ [ (-1„d“’‘i>o“g^(t)}  - o"’^DIJ®G^’(t)  + 0*'^1EJ“f^(t)}  - 0 (65) 

n i 


since  the  u”  are  orthogonal,  these  equations  uncouple,  and  we  ha/e 


Gn(")  + Vn<^>  - 


(66) 


lU 


where  determined  from  Equation  62^  is  found  to  be 

I / 1 


or 


Fn(T)  - I u’^Fdn  j I U®'^Du“dn 

Ojj  Oil 

Fn(T)  - j {Mv"  + Q<l>“}dn  / j U’^DU" 


(67) 

(68) 


o o 

The  modal  forces,  ^^(t),  may  be  computed  from  the  original,  physical, 

moments  and  forces, 

I 

11^  (n.T)«“(n)}dn 


(69) 


+($*^}dn 


Solving  the  set  of  ordinary  equations  (66)  for  the  scalar  time  dependent 
coefficients  Gj^(t)  and  substituting  Into  (63)  completes  the  solution  for 
the  forced  response.  The  method  of  Laplace  transforms  and  convolution  Is 
particularly  easily  to  apply.  If  the  initial  values  of  u,  u,  ^ and  ^ are 
all  zero,  then  the  initial  values  of  G^(t)  and  G^(t)  must  all  be  zero. 
Taking  transforms  of  Equation  (66)  then  yields 


Hence, 


and 


S g^(8)  + X^g^(8)  - f^(s) 
G„(t)  - f ^ (/^(T.p))F^(p)dp 


(70) 

(71) 


v*‘(nX 


)<l'(n,T)l  " 


T X, 

’ ^ (v^(T-6))f  '.-“(n)  + ♦“(!,)  }de 

J /J-  n J vkA  v<A 


to 


j + (♦")^)dn  (72)  J 

and  the  forced  response  Is  determined.  ° 

If  the  Initial  values  of  displacements  and  velocities  are  not  saro, 
the  free  response  may  be  determined  by  the  method  of  the  previous  section 
and  superposed,  or  the  initial  conditions  may  be  expanded  In  the  series 
given  as  Equation  (63),  0^(0)  and  G^(0)  determined,  and  Included  In  Che 
■Transformation,  Equation  (70).  A similar  result  has  been  given^. 
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IV.  Examples 


A.  Initial  Value  Problem. 

We  consider  a uniform,  simply  supported,  Timoshenko  beam,  subjected 

to  an  impulse  j,  per  unit  length  of  the  beam  at  time  t-0.  Hence,  at 
t^O*^,  the  beam  has  an  Initial  velocity,  given  by 


- I/pA 


Thus 


(73) 

k(x)  - l/pA,  f(x)  - g(x)  ■ h(x)  - 0 (74) 

The  appropriate  boundary  conditions  for  the  simply  supported  beam  are 


that 


♦ (o)  - ^(1)  - it  (o)  - (t)  • o 

•n  »n 


(75) 


To  find  the  eigenfunctions,  we  let 

(76) 

and  substitute  into  21  and  22  to  arrive  at  a characteristic  equation 

+ X[(l  + Y)p2  - 11  + Yp^  - o 


and  the  ratio  of  amplitudes: 


A/B  - (X  + p^)/p 


(77) 


(78) 


R>r  any  X,  there  are  four  values  of  p satisfying  equation  77,  which  we  de- 
note pj^  - ^^1»P2  “ P3  ■ ^^2*  P4  ■ "^’^2* 

- /{x(y  + 1)  - i)'^  + 4yX  }/2y 

- /{X(y  + 1)  + - 1)^  + 4yX  }/2y  (79) 

and  four  corresponding  ratios.  (A/b)^.Wlth  some  redeslgnatlon  of  terms,  we 
find 

• ■ bj^sln  rj^n  + b2Cos  rj^n  + b3Sln  r2n  + b^  cosr2n  (80) 


X-ri 


x-n 


f ■ -bj^  ( — j±-)cos  rj^n  + b2( — ^^)sln 

X-r?^  X-r,^ 

-bj  ( — ^ — ^)cos  r2n  + b^( ^)8ln  r2n 


(81) 
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The  boundary  conditions  for  the  simply  supported  beam  are  satisfied  by 
b2  ■ b^  » 0 and  either 


ril 

■ miT  , 

‘>3 

- 0 . 

m ■ 

1,2**- 

. . 

or 

^2^ 

■ nm  , 

*>1 

1 

O 

m ■ 

1,2*** 

• • •» 

Substitution  of  either  of  these  into  Equation  77  leads  to  the  same 
equation  for  X,  and  substitution  Into  Equations  80  and  81  leads  to  the 

same  mod a shapes;  hence  they  are  the  same  solution.  The  quantities  r^  may 
be  recognized  as  having  the  same  definition  as  the  wavenumber,  K,  IntrO" 

duced  previously.  Thus  we  see  that  motion  of  the  finite,  simply  supported 
Timoshenko  beam  can  occur  only  at  discrete  wavenumbers, 

K " r ■ mir/l  (82) 

For  each  m,  however,  there  are  two  eigenvalues,  X^  and  X^  2»  depicted  in 
Figure  2.  The  circles  In  that  figure  correspond  to  the  eigenvalues  for 
m • 1|2,3  and  4 In  a beam  for  which  1-10,  and  the  open  squares  denote 
the  eigenvalues  for  a beam  of  1 - 100  corresponding  to  m 1,  5,  15  and 
25.  From  this  figure  we  may  also  obseive  that  elementary  beam  theory  gives 
a good  approximation  to  the  first  5 frequencies  In  the  longer  beam,  but 
only  gives  a fair  approximation  of  the  first  in  the  shorter. 

Thus,  the  eigenfunctions  are 
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where 


2 r 2 2 

II  + (1  + Y)(”)  1 + [11  + (1  + y)(”) 


The  are  evaluated  with  the  negative  t gn,  end  evaluated  using 

the  positive.  From  these,  the  quantity 

I 

I u“Wd,  (85) 


■ay  be  evaluated,  with  the  result  ^ 

r nT  n r * ^>.1  1 a 

with  1«1  when  n-2m-l,  and  1-e  when  n-2n.  Substitution  into  Equations 
60  and  61,  and  the  result  Into  55  gives 


/ l-u . i co.:^  . 

(u/rl  B-l  I nr  ~ ^ 1 2 ") 

nim  L J 


m r 2 *T  — ^ (l-cos  ■ir)sln/x^  2 ’ 

^ r /mrv^  , c/X^,  ’ 

- “■  "" 

I < rnn/t  ?2 

sin  mvr)/l  J L ^ 

The  displacement,  u,  can  be  simplified  significantly. 


4Ir  7 
“ " TTw  2 


? sin  ntr  9 ct/iwxsto/x;;;^'  t> 


1,3,  •••-1-1  y j 

^ (mw/t)  y 


and  the  expression  for  the  rotation  is 


♦ ■-Mr  £ 


2 n I (s; 

j cos  mw  t ___? 


I 1(=1>  - V.tX.1-^ 


pAcl  i«i  /X  . /BE/  r 

IDfZ  'I.'  ' 


(mir/l) 


-}^  + 1 1 
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The  transverse  displacements  given  by  Equation  88  were  evaluated 


for  a bean  with  and  1-100.  The  peak  displacement  at  the  beam  center 

was  found  to  be  . 162.6  (^)  (90) 

occur Ing  at  t - 701.  For  a Bernoulli— Euler  beam  of  the  nnmr  geometry  and 
the  same  loading,  the  peak  response  was  found  to  be  within  IZ  of  this 

value,  and  occurred  at  t <■  681.  The  same  coiaputations  were  performed 
for  a shorter  beam  (I-IO) , with  the  results  that  a peak  displacement  of 

0)  = 1.756  (4Z£)  .... 

max  pAc  V.VJ-7 

was  found  for  the  Timoshenko  beam,  occuriMng  at  t -10.65,  and  8X  less 

for  the  Bemoul'il-Euler  beam,  occurring  at  t - 6.8. 

The  displacement  time  res^ionse  at  the  center  of  a beam  ri=10)  are 
compared  for  the  two  theories  in  Figure  3.  A more  dramatic  difference 
between  the  two  theories  occurs  in  the  predictions  for  the  bending 
moment.  Evaluating  the  first  derivative  of  Equation  89  at  the  beam 
center  (n  = 5)  leads  to  the  moment  time  history  shown  as  the  solid  line 
in  Figure  3,  while  the  second  derivative  of  the  series  for  the  displacements, 
i.e.,  the  moment,  as  determined  from  the  elementary  theory  does  not 
converge. 

B.  Traveling  Loads  on  the  Timoshenko  Beam. 

r 

Let  us  assume  that  a simply  supported  Tlmshenko  beam  la  aubj  acted 
to  an  external  load 

q(x,t)  - P«(x-vt)  (92) 

where  P Is  a force,  moving  at  apeed  v,  beginning  from  x-0  at  t-0.  We 
aasume  that  the  beam  dlaplacement  and  velocity  are  aero  at  t-0.  Here  P 
la  positive  In  the  -fa  direction  of  Figure  1,  In  terms  of  dlmensioolesa 
variables, introduced  through  Equation  19, 


19 


1 


and  m Is  taken  to  be  zero«  The  differential  equation  to  be  satisfied  is 


where 


^,nn  + ®“.n  + C“’°“,Tt+F-0 


■ vrt/c) 

U ■ < > and  F -<  ukA 


Ue  first  expand  the  applied  load  In  a series  of  the  eigenfunctions. 


F - X D F (T)a"(n) 
n“l 

and  look  for  the  response  In  terms  of  the  same  eigenfunctions, 

U - I G„(T)0"(n) 

The  response  Is  that  given  by  Equation  72. 

The  modal  forces  are  given  by  Equation  69,  or 


rP 


F„(t)  - 


{ — f 6(rn  - vrT/c)}sln  nir  dn 

JpUtcA ^ 

+ *n2}dn 


(HI)  _ xm,l 

,P  .miVTwri  .1 

F..(t)  “ t — r —r  — i'W  \ — * f'*’  ^ 

n'  ' yxA  1 c 2 I mrr  I 


where  the  choice  1*1  Is  again  used  for  n*2m-l,  and  1=2  Is  used  for  n*2m. 


The  complete  solution  Is 

- 2 j[(^) 


nn,T,  T «ir/t  V i sln222i TWTT— I 2 


cos  mffn/11 


•in[/x^(T-e)ide 


(101) 


Finally,  the  solution  for  displacement  may  be  evaluated  as 


{/ITT  sln«.lT  -Sf^^sln  ^Tl 
■ fl  l c 1 c *»* 


. . r ? rP  2 . n ® 

«(n,T)  - I I -Ter  T “Tv^rTTT? 1 — v; T 

^ (102) 


and  the  rotations  are 


^ P 2 . BIT  V rnir  V , n 

” ? :F=  Ai ^ =; 


wr/t 


(103) 


A series  representation  for  the  response  of  a Tlaoshenko  bean  to  a traveling 
load  has  been  previously  obtained  by  a different  method  . 

The  critical  speeds  for  the  traveling  load  occur  at 
the  vanishing  of  the  denominator  of  equation  102  , or 


V . ATTi- 


(104) 


c mir 

The  first  several  of  these  are  from  the  lower  branch.  For  l=10,  they  occur  at 


.163  (A^),.476  (iii),  .8125  (^),  1.147  (^),  etc. 
c » 2ir  3v  4v 

- .518  , .756  , .8621  , .913  , .940  , .957  (105) 

It  can  be  shown  that  these  values  assynptotlcally  approach  1.  The 
critical  frequencies  from  the  upper  branch  occur  at 

1-  1.211  (12),  1.657  (g).  2.187  (^),  2.753  (JJ),  etc. 

- 3.854  , 2.637  , 2.320  , 2.191  , etc.  (106) 

These  critical  values  asymptotically  approach  2.0.  or  in  general,  /-r.  For 

a longer  beam,  say  1=100,  the  critical  speeds  from  the  first  branch  occur  at 

I - .001969  (122),  .00782  (122),  .01739  (122),  .03042  (1^).  etc. 

^ It  2it  3ir  4v 


rlO^ 


rlOx 


r.l0. 


or  - - .0627  , .1244  , .1845  , .2421  , etc. 

c 


(107) 


and  asymptotically  approach  1 from  below.  For  the  critical  speeds  resulting 
from  upper  beanch  frequencies,  the  first  four  critical  speeds  are 


i 
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- = 1.002  (~),  1.010  (i^),  1.022  {—).  1.038  (~) , etc. 

C TT  bV  OJt 


(108) 


or  ^ = 31.89,  16.07,  10.84,  8.260,  etc. 

and  asymptotically  approach  2.0,  from  above,  or  in  general  /y. 

For  short  beams,  the  critical  frequencies  arc  densely  packed  in  two 
narrow  bands.  For  longer  beams,  they  are  somewhat  more  dispersed.  In  any 
case,  no  critical  speeds  occur  in  the  interval  1 £ v/c  £ /y. 

The  character  of  the  response  is  strongly  governed  by  the  ratio  of 
the  load  velocity  to  the  critical  velocity.  For  low  values  of  v (Figure  4) 
the  peak  displacement  occurs  ahead  of  the  point  of  load  application  at 
early  times.  For  a load  velocity  close  to  the  critical  value,  the  peak 
displacement  occurs  at  the  load  (Figure  5),  while  for  load  velocities  much 
higher  than  the  critical  value,  the  peak  displacements  are  seen  to  occur 
well  behind  the  load  (Figure  6).  In  each  case,  the  dot  denotes  the  posi- 
tion of  the  load  at  the  instant  the  displacements  are  given. 

The  shear  velocity  c,  is  also  seen  to  pl^y  an  important  role  in  the 
response.  In  Figures  4,  5 and  6,  the  location  of  a front  of  a wave  traveling 
at  speed  c is  denoted  by  the  short  vertical  line.  In  Figures  4 and  5,  low 
and  intermediate  velocities,  the  load  velocity  is  below  the  velocity  of  the 
shear  wave  and  no  significant  displacement  occurs  in  front  of  the  shear 
wave.  In  Figure  6,  the  load  velocity  is  much  in  excess  of  the  shear  velocity, 
and  the  peak  displacement  is  seen  to  occur  at  the  shear  front,  but  is  re- 
duced by  an  order  of  magnitude. 

It  is  of  interest  to  compare  these  results  to  the  response  of  a 
Bemoulli-Euler  beam  to  the  same  traveling  load,  so  that  the  influence  of 
rotatory  inertia  and  shear  deformations  on  the  motion  may  be  asaessed. 

The  differential  equation  to  be  solved  is  that  of  the  Bemoulli-Euler  beam,  or 
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+ pAu  • q ■ P6(x  - Vt) 


KIw 


xxxx 


(109) 


The  solution,  for  a simply  supported  beam  of  length  L,  initially  at  rest, 
is  found  to  be 


^ 2P  {sin  nitvt/L  - (nirv/Lu  )sin  u t}sin  mrx/L 
) fi-T- 5 0 n 

n-1  El  4 2_2  2/.  2 

^ n ^ - nx  V /L 


where  u ^ « EIn^ir^/(pAL^ 


) 


(110) 


(111) 


In  order  to  compare  the  results  of  the  elementary  theory  with  the 
results  from  the  Timoshenko  beam  theory,  we  make  the  following  substitutions 


T " t n “ */r,  f • L/r,  X/A  ■ r^,  • Ute/p,  as  before  and  set 

i . c2  - 

P K 

AT-  «„r/c  ■ /y 


2 2 


n n 


(112) 

(113) 


Thus 


u _ i£l  . i y fc  c t c ° 

“■  c2 

The  introduction  of  the  shear  coefficient  into  the  equations  is  totally 
■rtlflclal,  done  only  to  facilitate  comparison  with  the  Timoshenko  beam. 
The  critical  speeds  arise  at  the  vanishing  of  the  denominator,  and  are 


or 


^)  - .628  , 1.256  , 1.884  , 2.512  , etc.  (116) 

crlt 


These  velocities^  for  y=4  and  lal(\  increase  indefinitely,  rather  than 
reaching  an  asymptote,  as  was  found  in  the  case  of  the  Timoshenko  beam. 
For  tslOO,  the  critical  values  are  lOJt  of  those  given.  t 

Some  computed  results  are  given  in  Figures  7 and  8.  The  predictions 
of  the  elementary  theoiy.  Equation  214,  are  indistinguishable  (Figure  7) 
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0 .2  .4  *6  .8  1.0 

(a)  VELXITY  SUBSONIC,  BELOW  CRITICAL 


Figure  7 


Diaplaceaents  Doe  to  Traveling  Load,  Tiao^enko  or  Bemoulli- 
Euler  Bean  Theoiy,  1 > 100 
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A/2rP  - DISPLACEMENT  ^^A/2rP  - DISPLACEMENT 


from  the  results  of  the  Timoshenko  theory  for  the  long  beam  (L/r  = 100), 
even  for  speeds  near  to  the  first  critical.  In  the  shorter  beam  (L/r  = 

10),  however,  both  qualitative  and  quantitative  differences  become  apparent 
(Figure  8).  Since  the  elementary  theory  predicts  that  more  of  the  beam 
responds  at  early  times,  the  peak  displacement  is  significantly  under-^ 
estimated. 
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IV.  Summary 


I 


i 


i 

1 


An  orthogonality  relationship  for  the  Timoshenko  beam  has  been  shown 
to  make  possibl'i  the  solution  of  initial  value  problems  and  traveling 
load  problems  in  a manner  analagous  to  that  used  for  Bernoulli-Euler  beam 
problems.  Although  only  simply  supported  beams  were  considered  in  the 
examples  used,  the  orthogonality  condition  and  the  expansion  technique  is 
equally  applicable  to  any  of  the  four  possible  boundary  conditions. 

Through  an  example  problem  (the  simply  supported  beam  subjected  to  a uni- 
form impulse)  the  expression  for  the  moment  in  terms  of  Timoshenko  beam 
modes  was  found  to  converge,  while  the  moment  expressed  in  terms  of  the 
second  derivatives  of  the  displacements  of  the  Bernoulli-Euler  beam  did 
not. 

A traveling  load  problem  was  considered,  in  which  the  response  of  a 
simply  supported  Timoshenko  beam  to  a concentrated  transverse  force,  moving 
at  constant  speed,  was  determined  and  compared  to  analagous  results  for  the 
Bemoulli-Euler  beam.  For  the  Timoshenko  beam,  the  character  of  the  solu- 
tion was  found  to  depend  on  two  velocities  which  are  characteristic  of  the 
beam.  One  of  these  velocities  is  the  effective  speed  of  a shear  wave;  the 
other  velocity  is  that  velocity  of  the  traveling  load  which  would  lead  to 
unbounded  motion  in  a beam  of  infinite  length.  For  velocities  of  the  moving 
load  which  are  less  than  the  effective  shear  speed,  no  displacement  is  found 
ahead  of  the  shear  wave.  For  supersonic  loads,  i.e.,  loads  moving  at 
velocities  greater  than  the  effective  shear  wave  speed,  displacements  occur 
ahead  of  the  shear  wave  front,  but  the  maximum  displacement  occurs  at  the 
front.  For  very  low  values  of  the  speed  of  the  moving  load,  the  entire 
beam  is  involved  in  the  response,  and  the  response  resembles  motion  in  the 
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first  mode,  irith  the  peak  displacement  occurring  at,  or  near  the  beam 
center.  For  larger  values  of  the  traveling  load,  however,  the  peak  dis- 
placement  is  found  to  occur  tinder  the  load  unless,  • noted  above,  the 
load  is  supersonic. 

For  long,  thin  beams,  the  predictions  of  the  Bemoulli-Euler  theory 
do  not  differ  significantly  from  the  predictions  of  the  Timoshenko  beam 
theory.  For  shorter  beams  however,  significant  differences  are  found 
between  the  two  predictions.  The  predictions  of  the  elementary  theoiy 
do  not  show  the  previously  described  tendency  for  the  response  to  be  con- 
firmed to  the  portion  of  the  beam  processed  by  the  shear  wave,  for  such  a 
wave  is  not  'predicted  by  the  elementary  theory.  The  predictions  of  the 
elementary  theory  show  peak  displacements  at  or  near  the  beam  center  for 
velocities  of  the  traveling  load  where  the  more  complete  theory  predicts 
that  the  peak  displacement  occurs  under  the  load. 

A third  significant  difference  between  the  results  from  the  two 
theories  is  that  the  critical  velocities  of  the  Timoshenko  beam  Were  found 
to  have  a qualitatively  different  distribution  than  the  critical  velocities 
of  the  Bemoulli-Euler  beam. 


For  these  reasons,  it  appears  that  the  Timoshenko  beam  theory  should 
be  used  for  determining  the  response  to  traveling  loads  for  all  but  the 
most  slender  beams.  These  results  also  suggest  that  the  influence  of  shear 
deformations  and  rotatory  inertia  should  be  considered  in  determining  the 
response  of  more  complex  structural  elements,  such  as  rings,  plates,  and 
shells  to  moving  loads,  for  the  differences  between  the  predictions  of  the 
two  theories  for  beams  must  be  classed  as  qualitative,  rather  than  as  only 
qualitative.  Failure  to  consider  these  effects  may  lead  to  a failure  to 
model  important  physical  phemomena. 
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